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Measuring Expansion Velocities in Type II-P Supernovae 
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ABSTRACT 

We estimate photospheric velocities of Type II-P supernovae using model spectra 
created with SYNOW, and compare the results with those obtained by more conven- 
tional techniques, such as cross-correlation, or measuring the absorption minimum 
of P Cygni features. Based on a sample of 81 observed spectra of 5 SNe, we show 
that SYNOW provides velocities that are similar to ones obtained by more sophisticated 
NLTE modeling codes, but they can be derived in a less computation-intensive way. 
The estimated photopheric velocities (v mo d e i) are compared to ones measured from 
Doppler-shifts of the absorption minima of the H/3 and the Fe II A5169 features. 

Our results confirm that the Fell velocities (wj? e ) have tighter and more homoge- 
neous correlation with the estimated photospheric velocities than the ones measured 
from H/3, but both suffer from phase-dependent systematic deviations from those. The 
same is true for comparison with the cross-correlation velocities. We verify and im- 
prove the relations between VF e , vh/3 and v mo dei in order to provide useful formulae 
for interpolating/extrapolating the velocity curves of Type II-P SNe to phases not 
covered by observations. We also discuss the implications of our results for the dis- 
tance measurements of Type II-P SNe, and show that the application of the model 
velocities is preferred in the Expanding Photosphere Method. 
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1 INTRODUCTION 

Type II-P supernovae (SNe) are core-collapse (CC) induced 
explosions that eject massive hydrogen-rich envelope. These 
SNe show a nearly constant-luminosity plateau in their light 
curve. The plateau phase lasts for ~100 days, and finally 
ends up in a rapid decline of the luminosity as the expanding 
ejecta becomes fully transparent. 

CC SNe and their progenitors play important role in un- 
derstanding the evolution of the most massive stars. Exten- 
sive searches for SN progenitors showed that the progenitors 
of Type II- P SNe are stars with initial masses between ~ 8 
and 17 M jSmartt et alj|200gh . The lower value is roughly 
in agreement with the prediction of current theoretical evo- 
lutionary tracks, but the upper limit is somewhat lower than 
expected. The fate of the more massive (M ^ 20 Mq) stars is 
controversial. They are thought to be the input channel for 
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Type Ib/c SNe |Gaskell et allll986t ). b ut this has no t been 
confirmed directly by observations yet (Smartt 2009). 

Any observational study of the physics of SNe and their 
progenitors (and, obviously, all astronomical objects) de- 
pends heavily on the knowledge of their distances. Because 
distance is such a crucial parameter, strong efforts have been 
devoted to the development of distance measurement tech- 
niques that are applicable for Type II-P SNe. 

The Expanding Photosphere Method (EPM) 
(|Kirshner fe KwarJll974T ). a variant of the Baade-Wesselink 
method, uses the idea of comparing the angular and the 
physical diameter of the expanding SN ejecta. The input 
quantities are fluxes, temperatures and velocities observed 
at several epochs during the photospheric phase. Thus, 
EPM requires both photometric and spectroscopic monitor- 
ing of SNe throughout the plateau phase. A clear advantage 
of EPM is that it does not require external calibration via 
SNe with known distances. Another interesting property 
is that EPM is much less sensitive to uncertainties in 
the interstellar reddening and absor ption toward SNe. As 
lEastman. Schmidt fe Kirshnerl |l99^) (hereafter E96) have 
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shown, 1 mag uncertainty in Av results in only ~ 8 % error 
in the derived distance. However, the assumption that the 
SN atmosphere radiates as a blackbody diluted by electron 
scattering raised some concerns. These led to the develop- 
ment of the applications of full NLTE model atmospheres, 
like the PHOENIX code in the "Sp ectral- fitting Expan ding 



Atmosphere Method" ( SEAM) dBaron et al.1 |2004 ) 



the CMFGEN code by iDessart fc Hillierl l|2005al lbir The 
drawback of these approaches is that the building of tailored 
model atmospheres can be very time-consuming and needs 
much more computing power. Also, in order to compute 
reliable models, the input observed spectra need to have 
sufficiently high S/N and spectral resolution. 

The more recently developed S t andardized Candle 
Method (SCM) i|Hamuv fc Pintoll2002l ; iHamuv et al.ll2003l ) 
relies on the empirical correlation between the measured ex- 
pansi o n velocity and th e lumi nosity in optical (Nugen t et al 
20061; iPoznanski et all [2009) or in NIR (|Maguire et al 



2010a! ) bands in the middle of plateau phase, at +50 days 
after explosion. This method requires less extensive input 
data, but needs a larger sample of Type II SNe with indepen- 
dently known distances in order to calibrate the empirical 
correlation. Because this is basically a photometric method 
which compares apparent and absolute magnitudes, SCM is 
more sensitive to interstellar reddening than EPM/SEAM, 
as mentioned above. 

Beside photometry, both EPM/SEAM and SCM need 
information on the expansion velocity at the photosphere 
(vphot)- At first, EPM seems to be more challenging, because 
it requires multi-epoch observations, while SCM needs only 
the velocity at a certain epoch, at +50 days after explosion 
(1)50). However, direct measurement of V50 is possible only in 
the case of precise timing of the spectroscopic observation, 
which is rarely achievable. 

The correct determination of v p hot is not trivial. As 
SNe expand homologously (v ~ r, where v is the veloc- 
ity of a given layer and r is its distance from the center), 
in most cases it is difficult to derive a unique velocity from 
the observed spectral features. The most frequently followed 
approach relies on measuring the Doppler-shift of the ab- 
sorption minimum of certain spectral features (mostly Fell 
A5169 or H/3). Another possibility is the computation of the 
cross-correlation between the SN spectrum and a set of spec- 
tral templates, consisting of either observed or model spec- 
tra, with known velocities. The third method is building a 
full NLTE model (like PHOENIX or CMFGEN) for a given 
SN spectrum, and adopting the theoretical v p hot from the 
best-fitting model. 

The aim of this paper is to present a similar, but less 
computation-demanding approach to assign velocities to ob- 
served SN spectra. We apply the simple parametrized code 
SYN0W (|Fisherl Il999l ; lHatano et ail 1 19991 ) to model the ob- 
served spectra with an approximative, but self-consistent 
treatment of the formation of lines in the extended, homolo- 
gously expanding SN atmosphere. To illustrate the applica- 
bility of SYNOW, we construct and fit parametrized models to 
a sample of five well-observed, nearby Type II-P SNe that 
have a series of high S /N spectra publicly available. We also 
check and re-calibrate some empirical correlations between 
velocities derived from various methods. The description of 
the observational sample is given in fj2] In jj3] we first re- 
view the different velocity measurement techniques for SNe 



II-P, then we present the details of the application of SYNOW 
models ( i]3.4[) . The results are collected and discussed in i|H 
and while the implications of the results for the distance 
measurements are in ij6] We draw our conclusions in 33 



2 DATA 

Our sample contains 81 plateau phase-spectra of five ob- 
jects, SNe 1999em, 2004dj, 2004et, 2005cs and 2006bp, re- 
spectively. All five SNe are well-observed objects, they have 
been studied in detail, and show a wide variety in their phys- 
ical properties (see Table [TJ. 

SN 1 999em was discovered on 29th October 1999 by 
Q (|l999h in NGC 1637 at a very early phase. It is a 
very well-observed, well-studied objec t The explosion date 
was determined as 2451477.0 ± 2 JD dLeonard et all |2002| ; 



2006). We used the 



Hamuv et ail (|200ll ) 



IHamuv et al.1 l200ll; IDessart fc Hillierl 
spectra of iLeonard et al.l (|2002l ) and 
downloaded from the SUSPECTlj database. Our sample con- 
tains 22 spectra, covering the first 80 days of the plateu 
phase. 

SN 2004d,i was discovered on 31st July 2004 by Itagaki 
l|Nakano et alj|2004h in a young, massive cluster Sandage- 
96 of NGC 2403, ab out 1 month after explosion. We used 
the spectra taken by I Vinko et al.l |2006). Due to the lack of 
observed spectrophotometric standards, the flux-calibration 
of those spectra was inferior, but this is not a major concern 
when only velocities are to be determined. We included 12 
spectra taken between +47 and +100 days after explosion 
in our sample. 

SN 2004et was disco vered by Moretti 
l|Zwitter. Murani fc Morettil |2004| ) on 27th Sep tember 
2004 in NGC 6946. The spectra of ISahu et al.l (|2006[) 
(downloaded from SUSPECT) and iMaguire et al.1 (|2010bH 
were used, together with 6 previously unpublished early- 
phase spectra taken with the 1.88-m telescope at DDO (see 
the Appendix [XJ ■ The 22 spectra cover the period of +11 
to +104 days after explosion. 

SN 2005cs was discovered on 29th June 2005 bv lKloehrl 
(2005) in M51. Due to its very early discovery and proxim- 
ity, this object is very well-sampled and studied in detail. 
It is a underluminous, low-energy, Ni-poor SN that had a 
lo w-mass progenitor (see Table [T] for references) . 14 sp ectra 
of IPastorello et al.l (|2006l ) and IPastorello et all 12009) wort; 
used, which were obtained between days +3 and +61. 

SN 2006bp was d iscovered on 9th April 2006 by 
|Nakano fc Itagakill2006l ) in N GC 3953, also i n a ve ry early 
phase. We used 11 spectra of IQuimbv et al.1 l|2007l ). down- 
loaded from SUSPECT, covering the period between +5 and 
+72 days. 



3 THE ISSUE OF MEASURING EXPANSION 
VELOCITIES IN SNE 

During the photospheric phase the homologously expand- 
ing Type II-P SN ejecta consist of two parts: the outer, 
partly transparent atmosphere where the observable spec- 
tral features are formed, and an optically thick inner part, 
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Table 1. The physical properties of the supernovae used in this paper. 



SN 


to 


Distance 


cz 1 


E(B-V) 


M Ni 


Mprog 


References 2 




(JD- 245000) 


(Mpc) 


(km s _1 ) 


(mag) 


(10- 2 M ) 


(M©) 




SN 1999em 


1477.0 


7.5 - 12.5 


717 


0.10 


2.2 - 3.6 


<15 


1,2,3,4,5,6,7,8 


SN 2004dj 


3187.0 


3.2 - 3.6 


131 


0.07 


1.3 - 2.2 


12 - 20 


9,10,11,12,13,14,15 


SN 2004et 


3270.5 


4.7 - 6.0 


48 


0.41 


5.6 - 6.8 


9, 15 - 20 


16,17,18,19,20,21,22 


SN 2005cs 


3549.0 


7.1 - 8.9 


463 


0.05 


0.3 - 0.8 


6-13 


23,24,25,26,27,28,29 


SN 2006bp 


3835.1 


17.0 - 18.3 


987 


0.40 




12 - 15 


26,30 



NED, http://nedwww.ipac.caltech.edu/ 

2 References : (DlHamuv et al.1 120011) . (2) ILeonard et al.j (120021 ), (3) |Smartt et al.1 i2002h 



03]), (6) iBaron et al 
I d2004r j , (lOllKotak et al 



2053 ), (3)fSma_ 
2004h. (7) iDessart fc Hillierl toOoH ) 
j2005h . (ll)IWang et , al.l j2005n . ( 12)IChugai et al.. 
2i!06j), (15) IVinko et al.l J2009h . (161 [Li et al.l J2005I), 



(4) iLeonard et al 



2003), (5) Elmhamdi ct al. (200: 
2007 ). (9) lMaiz-ApeIIaniz et al 

200J ). (13) IZhang et al.l fcOOef), (14) fVinko et al.l 
(17) SahuetalJ (1200^^ (18) iMisra et al.1 i2007l), 
i 20M ), (21) iMaguire et al.l ll2010bl) , (22) ICrockett et al 
1 20051), (24) iPastorello et al l j2006l) , (2 5) iTakats 
Eldridge. Mattila fe Smartti | |2007T) 7(28) lUtrobin fc Chugail fcOOSl) . 



(8) [Ut robin 



1201 



Poznanski ct al 



Immler et al.1 J2007t) 



(2008) 
20091) 



Utrobin fc Chugail (120091) 

l|201ll) , (23) iMaund, Smartt fc Danzige' 
Vinkol (120061) . (26 ) IDessart et al.l 
(29) iPastorello et all 



(27) 
(30) 



which emits mo st of the continuum radiation as a "diluted" 
blackbody (e.g. E96). Because the inner part is mostly ion- 
ized, the major source of the total opacity is electron scat- 
tering. Thus, the instantaneous photosphere is located at 
the depth where the outgoing photons are last scattered, 
i.e. where the electron-sc attering optical depth is r e ~ 2/3 
(IDessart fc H illier 2005a], hereafter D05). This is presumed 
to occur in a thin, spherical shell at a certain radius r p hot, 
and the homologous expansion gives this a unique velocity 
Vphot — Vref ' Tphot /r re f, where r re / is the radius of an (arbi- 
trary) reference layer and v re f is its expansion velocity. As 
the ejecta expands and dilutes, the photosphere migrates 
inward, toward the inner layers of the ejecta that expand 
slower, thus, v p hot continuously declines in time. 

The issue of measuring v p hot comes from the fact that 
no measurable spectral feature is connected directly to v p hot- 
As previously mentioned in iJT] there exist several obser- 
vational/theoretical approaches to estimate v p hot from ob- 
served spectra. In the followings we briefly review these 
methods, summarizing their advantages/drawbacks. 

3.1 P Cygni lines 

The most widely used method is to measure the velocity 
represented by the minimum flux of the absorption compo- 
nent of an unblended P Cygni line profi le (sometimes re- 
ferred to as "maximum absorption" , v aos , IDessart fc Hillierl 
l2005bl ). This i s motivated by th e theory of P- Cygni line 
formation (e.g. K asen et al.l [2002). which predicts that the 
minimum flux in a P- Cygni line occurs at the velocity of 
the photosphere. Strictly speaking, this is only true for 
an optically thin line formed by pure scattering. In most 
cases the Fen A5169 feature is u sed for measuring v abs , 
for which IDessart fc Hilli er (2005b) pointed out that it can 
repre sent the true v P hot w ithin 5 - 10 % accuracy. How- 
ever, ILeonard et al. (2002) showed that using weaker, un- 
blended features one can get ~ 10 % lower velocities than 
from the Fen A5169 feature. This raises the question of 
whether the measurable lines were indeed optically thin and 
unbl ended, and which one repr esented the true v p hot- More- 
over, IDessart fc Hillierl l|2005bl ) also revealed that v a bs can 



either over- or underestimate v p hot, depending on various 
physical conditions, like the density gradient and the ex- 
citation/ionization conditions for the particular transition 
within the ejecta. 

While the Fen A5169 feature can be relatively eas- 
ily identified and measured in Type II-P SNe spectra ob- 
tained later than ~ 20 days after explosion (which is un- 
doubtedly a major advantage of this method), this is not so 
for earlier phases. In early-phase spectra only the Balmer- 
lines and maybe He I features (mostly the A5876 feature) 
can be identified. These features are less useful for measur- 
ing v p hot, because none of them are optically thin, and the 
He 1 A5876 feature can sometimes be blended with Nal D. 
Moreover, their line profile shapes often make the location 
of the minimum of absorption compo nent difficult to deter- 
mine (Des sart fc Hillier 2005bl. l2006l ). 

Also, in spectra having less signal-to-noise, the weaker 
Fe n features are often buried in the noise, making them in- 
appropriate for velocity measurement. In these cases sev- 
eral authors attempted to use the stronger feat ures, es- 
pecially Hp. Although IDessart fc Hillierl l2005bf ) pointed 
out that vh/3 is certainly less tightl y connected t o v v ho t 
than the weaker Fell A5169 feature, iNugent et all (2006) 
have shown that the ratio vhp/vfh is ~ 1.4 below 
vh/3 — 6000 km s , and it is linearly decreasing for 
higher velocities. Recently this correla tion was revised by 
IPoznanski. Nugent fc Filippenkc] (|2010l ) by revealing an en- 
tirely linear relation between vh/3 and VFe, but restricting 
the analysis only for velocities measured between +5 and 
+40 days after explosion. T he uncertainty of estim ating vfe 
from vh/3 is ~ 300 km s _1 (|Poznanski et al]|2010l ). It is em- 
phasized that these empirical relations are more-or-less valid 
between the velocities of two observable features, but both 
of th em may be systematically off from the true v p hot- 
deed. IDessart fc Hillierl (|2P05b) found that vh/3 can be either 
higher or lower than v p hot with an overall scattering of ±15 
%. For velocities below 10,000 km s _1 vnp is usually higher 
than v p hot, in accord with the empirical correlations with 
VFe, but above 10,000 km s" 1 it can be lower than v p hot- 
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-4000 -2000 2000 4000 -4000 -2000 2000 4000 

Velocity (km/s) Velocity (km/s) 

Figure 1. Model spectra of a Type II-P SN in early-phase (top left) and late- phase (top right), and the resulting CCFs (bottom panels) 
after cross-correlating the spectra with their identical models but having 2000 km s _1 higher v p h ot . The arrows mark the peak location of 
the CCFs, while the true relative velocity difference Avp^ot ls indicated by the dotted vertical line. The peak of the CCF underestimates 
Av p h ot in both cases. 



3.2 Cross-correlation 

Motivated by the uncertainties of measuring the Doppler- 
shift of a single, weak feature in a noisy spectrum, several au- 
thors proposed the usage of the cross-correlation technique, 
which is a powerful method for measuring Doppler-shifts 
of stellar spectra containing many narrow spectral features. 
This method predicts the Doppler-shift of the entire spec- 
trum by computing the cross-correlation function (CCF) 
between the observed spectrum and a template spectrum 
with a-priori known velocity vtempi- The resulting velocity 
is v cc — v re i + vtempi , where v re i is the velocity where the 
CCF reaches its maximum (i.e. the relative Doppler-shift 
between the observed spectrum and the template). 

However, the applicability of cross-correlation for SN 
spectra is not obvious, because of the physics of P Cygni 
line formation in SN atmospheres. If v p hot is higher, then 
the center of the absorption component gets blueshifted, 
while the center of the emission component stays close to 
zero velocity. Thus, cross-correlating two spectra with dif- 
ferent Vphot, the relative velocity will underestimate the true 
velocity difference between the two spectra. This is illus- 
trated in Fig. [1] where we determined the CCF between 
two Type II-P SN model spectra (computed with SYNOW, 
see below) that were identical except their v p hot which dif- 
fered by Av phot = 2000 km s' 1 . The cross-correlation was 
computed only in the 4500 - 5500 A interval to exclude 
the Ha feature. This was done for an early-phase spectrum 
containing only H and He features (Fig. [1] left panel) and 
a later-phase spectrum with more developed metallic lines 
(Fig. [T] right panel). It is visible in the bottom panels that 
the maximum of the CCF (i.e. v re i) is shifted toward lower 



velocities with respect to the true Av p hot (indicated by a 
dotted vertical line in the CCF plot) for both spectra. The 
systematic offset of v re \ from Av p hot is ~ 200 km s _1 for 
the later-phase spectrum that contains many narrower ab- 
sorption features, and it reaches ~ 700 km s _1 for the early- 
phase spectrum that are dominated by broad Balmer lines 
with stronger emission component. The se simple tests are in 
ver y good agreeme nt with the results of lHamuv et al.l (1200 ll ) 
and lHamu vf l|2002l ), who applied the cross-correlation tech- 



nique using model spectra of E96 with a-priori known v p hot 
as templates. They pointed out that v re i usually underesti- 
mates Avphot so that Avphot/v T ei ~ 1.18, with the scatter 
of ~ 900 km s" 1 . 

Another varian t of th e CCF method was proposed 
by IPoznanski et ail (20091) . and it was als o app ied by 
ID'Andrea et al.l l|2010l) and IPoznanski et all l|2010l ). They 
took their template spectra from the library of the Su - 
perNova IDentification code SNIE0 (|Blondin fc Tonry|[2007l ) 
that contains high S/N observed spectra of many SNe, but 
used only those templates that showed well- developed Fell 
A5169 feature. The reason for their choice was that they 
wanted to estimate VFe at +50 days after explosion (vso), 
which is an input parameter in SCM. This template selection 
caused the well-known issue of template mi smatch that fur- 
ther b iases the cross-correlation results. As ID'Andrea et all 
(|2010h concluded, this template mismatch can result in sig- 
nificant underestimate of VFe by ~ 1500 km s _1 for spec- 
tra obtained at t < 20 days, and is also present in later- 
phase specta, although being less pronounced, ~ 300 — 400 
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km s _1 . To over c ome this problem, in a subsequent paper 
IPoznanski et alj |2010l ) suggested the application of their 
VFe vs. VHp relation to estimate VFe by measuring vh/3 for 
these early-phase spectra and then propagate the resulting 
velocity to day +50. Nevertheless, the two groups presented 
velocities that were different by 200 to 1000 km s _1 for the 
same set of SNe spectra from the SDSS-II survey. This un- 
derlines that although cross-correlation seems to be an easy 
and robust method that gives reasonable velocities even for 
low S/N spectra, the results may be heavily biased, espe- 
cially at early phases, when the SN spectra contain only a 
few, broad features. 



3.3 Tailored modeling 

Tailored modeling of the whole observed spectrum via 
NLTE models were invoked e.g. by baron et alj (I2004T ) us- 
ing the code PHOEN IX, and iDessart fc Hillierj |200rj ') and 
iDessart et al. (2008) with the code CMFGEN. Here v p hot is 
determined implicitly by requiring an overall fit of the en- 
tire observed spectrum by a synthetic spectrum from a full 
NLTE radiation-hydrodynamics model. In these models the 
location of the photosphere is usually defined as the layer 
where t he ele ctron scattering optical depth reaches unity or 
2/3 (e.g.|D(l). Vphot is then simply determined by the radius 
r p hot and the law of homologous expansion (see above) . Al- 
though this is probably the best self-consistent method for 
obtaining SN velocities, building full NLTE models for mul- 
tiple epochs requires a lot of computing power. Thus, its us- 
age for a larger sample of SNe would be very time-consuming 
and impractical. Also, in some cases a good global fit to 
the whole spectrum may not be equally good for individual 
spectral lines, where lots of physical details play important 
role in the formation of the given features. This may lead 
to some uncertainties in the velocities given by the models, 
which will be illustrated in §2] 



3.4 Using SYNOW 

The success of tailored spectrum models to estimate v p hot 
suggests a similar, but much more simplified approach: the 
usage of parametrized spectrum models that do not contain 
the computation-intensive details of NLTE level populations 
or exact radiative transfer calculations, but preserve the ba- 
sic physical assumptions of an expanding SN atmosphere, 
and able to reproduce the formation of P Cygni lines in a 
simplified mann er. Such models can be computed either with 
the SYNOW code (Branch et al. 2002) , or the more recent ly de- 
veloped SYNAPPS code l|Thomas. Nugent, fc Mezal201ll ) that 
has a parameter-optimizer routine built-in. In this paper we 
apply SYNOW to calculate model spectra. Because SYNOW does 
only spectrum synthesis and has no fitting capabilities, we 
used self-developed UNIX shell scripts to fine-tune the pa- 
rameters until a satisfactory fit to the observed spectrum is 
achieved. 

The basic assumptions of SYNOW are the followings: i) 
the SN ejecta expand homologously; ii) the photosphere ra- 
diates as a blackbody; Hi) spectral lines are formed entirely 
above the photosphere; iv) the line formation is due to pure 
resonant scattering. Level populations are treated in LTE, 
and the radiative transfer equ ation is solved in th e Sobolev 
approximation (see also in e.g. Ka sen et al.l 12002). 



When running SYNOW, several parameters must be set. 
These are the temperature of the blackbody radiation (Tbb) 
emitted by the photosphere, the expansion velocity at the 
photosphere (v p hot), the chemical composition of the ejecta 
and the optical depth of a reference line (i~ re /) of each com- 
pound. For each atom/ion the optical depths for the rest 
of the lines are calculated assuming Boltzmann excitation 
governed by the excitation temperature T exc . The location 
of the line-forming region in the atmosphere can be tuned 
for each compound by setting the velocities of the lower and 
upper boundary layers, v m i n and v ma x- The optical depth as 
a function of velocity (i.e. radius) can be modeled either as a 
power-law, or an exponential function. We assumed power- 
law atmospheres, and adjusted the power-law exponent n to 
reach optimal fitting. 

After setting the inital values by hand, several models 
in a wide range of v p hot, n, and r re f were created for a 
pre-selected set of ions. In order to reduce the number of 
free parameters, we initially set Ttt (which has very little 
effect on the line shapes) to represent the continuum of the 
fitted spectrum and kept it fixed during the optimization. 
Moreover, we applied a single power-law exponent n for all 
atoms/ions. We also assumed that all spectral features are 
photospheric, thus, fixing v m in well below the photosphere 
and v m ax at ~ 40000 km s _1 . 

The best-fitting model was then chosen via \ 2 ~ 
minimization, and the fitting process was iterated for a few 
times, each time resampling the parameter grid in the vicin- 
ity of the minimum of the \ 2 function found in the previous 
iteration cycle. This way we determined the parameters and 
the chemical compositions that best describe the observed 
spectra. 

Then, to further refine the estimated photospheric ve- 
locity, we fine-tuned only v p hot of the best-fitting model, 
and calculated the \ 2 function only in the vicinity of cer- 
tain lines instead of the whole spectrum. This may reduce 
the systematic under- or overestimate of v p hot produced by 
false positive fitting to the observed spectrum outside the 
range of the considered spectral features. 

Motivated by the results of IDessart fc Hillierj (|2005bh 
(see i]3.1[) . we have chosen the Fell A5169 feature for this 
fine-tuning process. When this feature was not present in 
the observed spectrum (i.e. the early-phase spectra, before 
~15 days) we used H/3 (see i)3.1[) instead. Hereafter we de- 
note the Vphot parameter of the best-fitting SYNOW model as 
Vmodei ■ Errors of v mode i were estimated by choosing the 90 % 
confidence interval around the minimum of the \ 2 function. 

Fig. [2] shows two examples for an early- and a later- 
phase spectrum of SN 1999em together with the best-fitting 
model. The right panel zooms in on the region of H/3 and Fe II 
A5169. Note that although the final fitting was restricted to 
the proximity of these lines, the best model describes the 
entire observed spectrum (except Ha) very well. 

This velocity measurement method have multiple 
sources of error. One of them may be the systematic bias 
due the approximations in the model (LTE, power-law at- 
mosphere, simple source function, etc.). However, the com- 
parison of our results with those from full NLTE CMFGEN 
models (Q show no systematic bias in the case of SNe 
1999em and 2005cs. The agreement between the velocities 
from these two very different modeling codes are within ±10 
percent. For SN 2006bp the differences are higher, but it will 
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Figure 2. The observed (black solid line) and best-fitting model (red dashed line) spectra of SN 1999em on days +9 (top) and +41 
(bottom). The right panels enlarge the wavelength regions of the fitted lines: H/3 (top) and Fell A5169 (bottom). 



be shown below that for this SN the CMFGEN models do not 
describe well the spectral features we use, contrary to the 
SYN0W models (g3J. 

Another source of error may be the correlation between 
the parameters. In Fig. [3] we present contour plots of the 
X 2 hyperspace around its minimum, as a function of v mo d e i 
and several other parameters that can affect the shape of 
the fitted Fell A5169 feature. The thick black contour curve 
corresponds to 50 % higher \ 2 than the minimum value. It 
is visible that correlation is indeed present (i.e. the contours 
are distorted) between v mo dei and the power-law exponent 
n or the optical depth tf r - The correlation is much less 
between v mo d e i and r re / of Tin and Mgl, whose features 
may blend with Fe II A5169. However, even for the correlated 
parameters, selecting n or TFe very far from their optimum 
value can alter v mo j e i only by a few hundred km s _1 . Thus, 
we conclude that uncertainties in finding the minimum of \ 2 
do not cause errors in v mo d e i that significantly exceed the 
uncertainty due to the spectral resolution of the observed 
spectra (which is usually between 200 - 300 km s _1 ). 

A possible source of uncertainty may be that during the 
final fitting the wavelength interval around the used spectral 
feature is chosen somewhat subjectively. However, our tests 
showed that changing the limits reasonably has negligible 
effect on the final velocities. 

It is emphasized that although the final fitting is re- 
stricted to a vicinity of a well-defined spectral line, this 
method is certainly more reliable than the measurement of 
only the location of the minimum of the same feature. As it 
was discussed above, the minimum can be significantly and 
systematically altered by signal-to-noise, spectral resolution, 
blending, etc. The fitting of a model spectrum to the entire 
feature is expected to overcome these difficulties, provided 
the underlying model is not too far from reality. 

4 COMPARING THE RESULTS FROM 
DIFFERENT METHODS 

Using SYN0W as described above, we determined the best- 
fitting parameters of all SNe spectra from Sec. [2] The result- 



ing model velocities are collected in Table [BT1 in Appendix 

B. The best-fitting SYNOW parameters, such as T re / for each 
atom/ion, the power-law exponent n and v mo dei together 
with the chosen T p hot, can be found in Table lCll in Appendix 

C. In Table iBll we also list the vf b and vh/3 velocities. For 
SNe 1999em, 2005cs and 2006bp, we col lected the photo- 
spheri c vel ocities from CMFGEN models of iDessart fe Hillierl 
(2006) and lDessart et all l|2008h . These are included in Ta- 
ble IBll as v n ite ■ Velocities from the cross-correlation tech- 
nique (Sect. 13. 2[) were obtained using two sets of template 
spectra. The first set contained the 22 observed spectra of 
SN 1999em (set #1), while the second set was based on the 
CMFGEN models mentioned above (set #2). The velocities of 
the template spectra were VFe for set #1 and v n ite for set 
#2. We cross-correlated all the observed spectra with the 
two sets separately on the wavelength range of 4500 - 5500 
A, and the resulting velocities are also in Table [BT1 as v cc . 

Fig. [4] shows Vmodei against phase for all studied SNe 
(top left panel), and the ratio of v mo d e i to all the other 
velocities. The calculated velocities all show the expected 
decline with phase as the photosphere moves deeper and 
deeper within the ejectra, toward slower expanding layers. 

Similar plots containing the ratio v mo d e i/v cc and 
Vabs/vcc as functions of phase, are presented in Fig. [5] 

In the followings we provide some details of deriving 
these velocities for each object and discuss some object- 
specific differences between them. 



4.1 SN 1999em 

When determining Vmodei with SYNOW, H/3 was fitted for the 
first 6 spectra, then the Fell A5169 feature was used for the 
remaining 16 spectra. The resulting velocities are between 
11000 and 1800 km s" 1 . As seen in the bottom right panel 
of Fig. 13] v mo dei and Vnp are about the same for the early 
phases (before the appearance of the Fell lines), while later 
vnp tends to be higher than v mo dei. Also, between day +15 
and day +40, v mode i is a slightly higher than vf s (Fig. 2] 
bottom left panel). After day +40 v mo d e i drops below vfb 
and their ratio increases toward later phases. 
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Figure 3. Contour plot of the x 2 function around its minimum, as a function of v mo[ i e i and power-law exponent n (a), r re f of Fen (b), 
Till (c) and Mgl (d). The thick black contour curve corresponds to 50 % higher x 2 than the minimum value. The parameters plotted in 
(a) and (b) are definitely correlated, while the correlation is much less between the parameters in the (c) and (d) panels. 



The velocities from CMFGEN models of (|Dessart fc Hillierl 

2006) (Fig. [4] top right panel) agree with v mo dei- The cross- 
correlation with set #2 (FigO bottom panels) gave similar 
results for the first few points, but overestimate v m odei be- 
tween days +22 and +80. They mostly fall between imp and 
VFe, which is expected, since we cross-correlated the range 
of 4500 - 5500 A, where these features appear. 



4.2 SN 2004dj 

The SYN0W model velocities of the 11 spectra that cover the 
second half of the plateau phase are between ~ 3400 and 
1700 km s _1 . These are similar to those of SN 1999em at 
the same phase. Both vfb and vnp are higher than Vmodei 
at all epochs, especially the latter with a factor of about 1.8 
(Fig.©. 

No CMFGEN model was available for SN 2004dj. Cross- 
correlation with both template sets gave very similar re- 
sults. They are only slightly higher than both v mo d E i and 
v Fe (Fig[5j. 

4.3 SN 2004et 

For the first 6 spectra the SYNOW model was optimized for 
Hf3, then for the Fell A5169 feature. The resulting model 
velocities are between 9700 and 1800 km s" 1 (Fig. gj). The 



VFe values are similar to v mo del, but their ratio shows slight 
phase dependence, similar to the other SNe studied here. 
On the contrary, the values of Vnp are very different from 
Vmodei ■ At early phases they are close to v mo dei (except for 
the first point), but later the vup to v mo dei ratio strongly 
increases, becoming as high as 2.5. 

Again, there is no CMFGEN model available for this SN. 
Cross-correlation with set #1 resulted in velocities similar 
to vnp at early phases and to VFe later. With set #2, cross- 
correlation gave similar results at early phases, but later it 
produced systematically higher velocities. This underlines 
the importance of selecting proper template spectra and 
template velocities when applying the cross-correlation tech- 
nique. 

4.4 SN 2005cs 

We used the H/3 line for fitting the first 3 spectra with SYNOW. 
The velocities of this SN are very low: they are in the range of 
7100 - 1100 km s _1 and decrease quickly. The velocities from 
absorption minima are very close to v mo dei for both H/3 and 
Fell A5169. The VFe values follow the tendency similar to 
the previous objects: they are somewhat lower than v mo dei 
at the early phases, but get higher after about day +30. 
The Vup values are much closer to v mo d e i than for the other 
SNe, and the VHp/Vmodei ratio stays about the same for all 
epochs (Fig. 0}. The velocities of the CMFGEN models for 
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Figure 4. Top left panel: model velocities (v mo d e i) from SYNOW as functions of phase. Top right panel: the ratio of SYN0W and CMFGEN model 
velocities (v mo d e i / v n i te ) against phase. Bottom panels: The phase dependence of vp e /v mo d e i (bottom left) and v H p/v rno d e i (bottom 
right). Different symbols code the followings: filled circles - SN 1999em; open circles - SN 2006bp; filled triangles - SN 2005cs; open 
triangles - SN 2004dj; asterisks - SN 2004et. 



SN 2005cs are the same as the v mo d e i values for all epochs, 
except for day +9. Cross-correlation with both template sets 
resulted in velocities close to vp e . 

4.5 SN 2006bp 

The results for SN 2006bp are controversial. Applying SYNOW, 
the H/3 line was fitted for the first 4 of the 11 observed 
spectra, while Fen A5169 was used for the rest. The model 
velocities are between 12000 and 3800 km s" 1 . Both VFe and 
Vh/3 follow the tendency shown by other SNe (Fig. 2J. 

On the contrary, unlike in the previous two cases, the 
velocities from the CMFGEN models differ significantly from 
our v mo dei values. At early epochs this difference is much 
lower (~ 500 — 700 km s _1 ) being close to zero at day +15. 
After day +15 it gets higher reaching ~ 1200 km s" 1 on day 
+32. At later phases the difference decreases somewhat, but 
stays being significant. 

Cross-correlating the same spectra with the CMFGEN 
models using the wavelength range of 4500 -5500 A resulted 
in velocities that are very close to v mo d e i (except for day +9) . 
Using the #1 template set, the results agree well with dfc 
or vh/3 at early epochs. 

To examine the obvious controversy between the veloc- 
ity of the CMFGEN models and all the others, we plotted the 
observed spectra and the best-fitting CMFGEN model on day 
+32 (when the differences are the highest) in Fig. [6] Zoom- 
ing in on the range of 4500 — 5500 A clearly shows that the 



model by Dessart et al.l (|2008l ) does not fit these spectral 
features well, leading to an underestimate of the velocity. 
Thus, we suspect that the velocity differences we found are 
probably due to the inferior fitting of the CMFGEN models to 
the SN 2006bp spectra. 



5 DISCUSSION 

As shown in the previous sections, the photospheric veloci- 
ties of four SNe in our sample evolved similarly. SNe 1999em, 
2004et and 2006bp had high velocities at early phases and 
they decreased quickly, although their decline slopes were 
different. SN 2004dj probably showed similar evolution, but 
the lack of the early-phase data prevents a more detailed 
comparison. On the contrary, SN 2005cs was a very differ- 
ent, low-energy SN II-P as discussed in detail in previous 
studies. It had lower early velocities and the velocity curve 
decreased much faster than for all the other SNe. 

As expected, the different velocity measurement meth- 
ods we applied provided somewhat different results. As seen 
in Fig[5] the velocities obtained from cross-correlation are 
usually closer to v a b s than to v mo d e i ■ This is understandable, 
given that the cross-correlation method is most sensitive to 
the shapes and positions of the spectral features that may 
be biased toward lower, or higher velocities. The v mo d e i/vcc 
ratio (Fig[5]top left and bottom left panels) shows the same 
trend (but plotted upside down) as the VFe/vmodei ratio in 
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Figure 5. Top panels: Ratio of f mo( jei an d ^abs t° "cc from cross-correlating with the observed spectra of SN 1999em (set #1, see text) 
as functions of phase. Bottom panels: the same as above, but with respect to the template set containing CMFGEN models (set #2). 
The symbols code the same SNe as in Fig|4] 
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Figure 6. Plot of the CMFGEN (red dotted line) and the SYNOW models (turquoise dashed line) and the observed spectra (black 
continuous line) of SN 2006bp on day + 32 (see text). 



FigH] (bottom left panel), i.e. v mo d e i is higher between day 
10 and 50, but becoming smaller than vf b or v cc . On the 
other hand, no such systematic trend can be identified be- 
tween v m odei and v n ite (Fig[4] top right panel). These bench- 
marks suggest that the model velocities, either from SYNOW or 
CMFGEN are consistent, and they show phase-dependent off- 
sets from the absorption minima, or cross-correlation veloci- 



ties. The increasing systematic offset is particulary strong for 
v Hfi (FigH] bottom right panel). Thus, the traditional, sim- 
ple measurement methods seem to underestimate the true 
photospheric velocities before day 50, but increasingly over- 
estimate them toward later epochs. This should be kept in 
mind when the true photospheric velocities are needed, e.g. 
in the application of EPM. 
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Figure 7. The vpe/^model (left panel) and vh8 l v model (right panel) ratio against the model velocities. The symbols have the same 
meaning as in Fig|4] 



In order to do further testing, we plotted the ratio 
of v Fe /v mode i and Vh g /vmodel as a functio n of v mo dei for 
all SNe, following IPessart fc Hiliien (|2005bl ) (Fig. 0. The 
VFe/ ' Vmodel ratio shows the same trend for all objects: at high 
velocities (i.e. early phases) the ratio is somewhat lower than 
1, then it reaches unity around ~ 4000 km s _1 , and below 
that it keeps rising, reaching ~ 1.6 by the end of the plateau. 
The Vh g / Vmodel ratio is more complicated. At high v mo dei 
values it is around 1, but becomes higher than unity around 
Vmodei ~ 7000 km s _1 . Below that the slope of the rising 
changes from object to object. In the case of SNe 2005cs 
and 1999em this ratio stays under ~ 1.4, while for the other 
three SNe it becomes much higher. For SN 2006bp there are 
no spectra below v mo dei = 3850 km s _1 , but above that its 
evolution seems to be similar to that of SN 2004et. 



A similar plot was published by lPessart"fc Hillicr 



(2005b) based on their set of CMFGEN model spectra (see 
their Fig. 14). The only slight difference is that they plot- 
ted the ratio of the velocity measured from the absorption 
minima of the model spectra to the input velocity of the 
code, as a function of the input velocities. Although they 
did not have data below ~ 4000 km s _1 , and we do not 
have data above ~ 12050 km s _1 , between these limits their 
plotted values are mostly similar to ours. In their Fig. 14 the 
Fell A5169 velocities are lower than that of the model for 
high velocities, and their ratio reaches 1 between 5000 and 
4000 km s" 1 , just like our data. The situation is somewhat 
different for H/3. At high velocities the two results are consis- 
tent: ab ove ~ 11000 km s _1 the data bv lDessart fc Hillierl 
(2005b), as well as ours, are around 1. However, their ve- 
locity ratio exceeds 1 at ~ 8000 km s^" 1 and has a high- 
est value of 1.15 for H/3. It is much lower than our results 
in Fig0 In the case of SN 2004et, our velocity ratio goes 
as high as 2.5. I t must be noted, howeve r, that the model 
spectra used bv lDessart fc Hillierl ( 2005bl) were tailored to 
represent SNe 1987A and 1999em (|d05| ). The latter object 
is also in our sample, and our vh8 to v m o dei ratio for that 
particular SN is similar to the results of Pess art fc Hillierl 



(|2005bl). Thus, it is probable that the lower v H 8/ v modei ratio 
of Pessart fc Hil licr ( 20053) is due to the limited parameter 
range of their CMFGEN models used to create their plot. 



Recentlv lRov et al.l (|201ll ) published a study of velocity 
measurement for the Type II-P SN 2008gz. They applied a 
similar technique of using SYNDW to fit Fen features around 
5000 A. They also estimated the velocity from the absorp- 
tion minima of these lines. They got 4200 ±400 km s" 1 and 
4000 ± 300 km s _1 for v mo dei and VFe, respectively, from a 
+87d spectrum. This result is consistent with our findings 
plotted in Fig. [7| vf e is practically equal to v mo dei around 
4000 km s" 1 . 



5.1 Velocity-velocity relations 

Using the synthetic spectra of |E96| and |P05l . iJones et al.l 
(2009) also examined the relation between vh8 and Vmodel- 
They found that their ratio can be described as 



VH8 



E 

3=0 



a 3 V H8 



(1) 



where the values of aj are given in Table [2] In Fig. 8(a) 



plotted our data toget her w ith these polynomials. The poly- 
nomials based on the models overestimate our v mo dei 
values (rms a = 0.412), while those from the E96 models 
provide much better fit for all SNe except SN 2004et (rms 
a = 0.301, but a = 0.178 without the d ata of SN 2004et). 
We fitted Eq.JTJ to our data (Fig. 



black curve) 



The resulting aj coefficients are in Table El Our fit resulted 
in a much lower rms scatter, a = 0.276. Repeating the fitting 
while omitting the data of SN 2004et, the result became very 
similar to that from the E96 models. 

Since VFe is thought to be a better representative of 
the velocity at the photosphere than Vbb, it is expected 
that v m odei can be predicted with better accuracy by mea- 
suring VFe- Indeed, Fig. [7]suggests that the VFe/vmadel ratio 
is almost the same from SN to SN, unlike the vnp/vmodei 
ratio that can be quite different for different SNe. Thus, we 
repeated the fitting of Eq.{T} using VFe instead of Vhb (Fig. 
8(b) |. We found the rms scatter of a = 0.111, which is much 



lower than in the previous cases. The aj coefficients of this 
fitting are also included in Table [2] 
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Figure 8. (a) The plot of t> m odei ag ainst vjjp. The lines represent the polynomials calculated by I Jones et al, | (1200911 (see text) based on 
the models of lE96l (red dotted line) . [Doil (blue dashed line). The result from our fitting on all SNe is plotted as a black dotted line line 
while the fitting omitting SN 2004et is shown by a black solid line, (b) The same fitting as (a) but for vp e using all SNe. The simbols 
have the same meaning as in Fig|4] The polynomial coefficients are in Tabic [2] 



The tight relation between vfc and v mo dei in Fig. 8(b) 
suggests a possibility to estimate v mo d e i from the measured 
VFe values. However, it is emphasized that SN-specific dif- 
ferences in the expansion velocities may exist, thus, model 
building for a particular SN, whenever possible, should al- 
way s be preferred . 

INugent et alJ l|2006l ) found that VFe evolves as 



VFe(t)/v F e(50d) = (i/50) C 



(2) 



where c = —0.464 ± 0.017. After repeating the fitting of 
Eq.Q to our data, we found the exponent to be c = 
-0.663 ± 0.01. Then, since the data of SN 2005cs are very 
different from the rest of the sample, we omitted the ve- 
locities of SN 2005cs and repeated the fitting. This re- 
sulted in c = -0.546 ± 0.01 (Fig. |9(a)[ ). These two expo- 
nents marginally di ffer (at ~ la) from the value given by 
iNugent et all (|2006h . A possible source of this difference (be- 
side the different velocity measurement techniques applied) 
may be that our sample covers the phases be t ween +13 and 
+ 104 days, while the data by INugent et al.1 (|2006T l are be- 
tween +9 and +75 days. 

We also examined how the SYN0W model velocities evolve 
in time. Combining Eq.JT]) and Eq.©, the following relation 
has be derived (again, excluding SN 2005cs from the sam- 
pie): 



(t) /v m odel (50d) 



(t/50)- 



E 



(3) 



where b = 0.467 ± 0.15, 6i = 0.327 ± 0.23 and b 2 = 0.174 ± 



0.11. The rms scatter is a = 0.148 (Fig. [9(b) I. 

As was mention ed in £13.11 using SDSS data 
IPoznanski et al.l (|20101 ) examined the correlation between 
velocities measured from the absorption minima of H/3 and 
Fen A5169 lines (see Fig. [10]). They found that there is a 
linear relation given by VFe(50d) — a ■ ^H,9(50d), where 
a = 0.84 ± 0.05. Using our sample we repeated their fit- 
ting. First, we used all epochs where both vnp and vfh were 
measured. The slope of the fitted line was a = 0.791 ± 0.012 
(a — 0.146). Then, we kept only the velocities obtained be- 



Table 2. Polynomial coefficients for the vjja and vp e to f mo< j e ! 
ratio (Eq[T}. 







1 



ref. 



aj (Hp) (E96) 1.775 -1.435e-4 6.523e-9 0.30 (1) 

aj (Hp) (D05) 1.014 4.764e-6 -7.015e-10 0.41 (1) 

a 3 (Hp) 1.528 -1.551e-5 -3.462e-9 0.27 (2) 

a 3 (Hp) w/o 04et 1.578 -8.573e-5 3.017e-9 0.17 (2) 

aj (Fen A5169) 1.641 -2.297e-4 1.751e-8 0.11 (2) 



(1) Ijones et al] l|2009l 'l 

(2) this paper 
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Figure 10. The relation between the measured vp e and Vjjp 
values. The empty circles show the data taken after day 40, while 
the filled circles refer to data taken before day 40. The line shows 
the fit to the filled circles only. 



fore day 40 (similar to IPoznanski et aT1l2010l ). This resulted 
in a — 0. 823 + 0.015 (a = 0.102), which is basically the same 
as that of IPoznanski et all (|2010l). Thus, ou r study fully con- 
firms the results by Poznanski et al.l j2010r i. but extends the 
validity of the relation toward later phases. 
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6 IMPLICATIONS FOR DISTANCE 
MEASUREMENTS 

Improving the accuracy of the velocity measurements has an 
important aspect in measuring extragalactic distances with 
SCM or EPM (cf. §1). Both SCM and EPM need velocities, 
thus, the results of this paper can be significant for both 
techniques. 

SCM was calibrated using vfb on day +50. However, 
in many cases getting a spectrum at or around day +50 
is not possible. In these cases Eq.© can be used to esti- 
mate i>Fe(50d). We have improved the exponent in Eq.([5]) 
as —0.546 ± O.Of, based on more data obtained on wider 
range of phase than previously. The difference betw een our 
result and the previous curve l|Nugent et al.l 120061 ) is the 
highest around day +20 (Fig[9}. The new curve may result 
in better constrained iiFe(50d) when only early- phase spec- 
tra obtained around day +20 are available. 

However, there are several drawbacks of SCM. For ex- 
ample, the uncertainty in the moment of explosion, i.e. in de- 
termining the phase of a particular spectrum can lead to sig- 
nificant error in the distance determination. Moreover, as the 
example of SN 2005cs shows in FigO some Type II-P SNe 
can deviate significantly from the average, especially during 
early phases. Thus, one should be careful when such kind 
of interpolation or extrapolation is to be applied. The ex- 
ample of SN 2005cs suggests that multi-epoch spectroscopic 
observations should always be preferred against single-epoch 
spectra when distance determination is the aim. 

The case of EPM is different. Since this method does 
not require calibration, but needs multi-epoch data, devia- 
tions in the measured velocities have higher impact. To show 
this, we calculated t he EPM-dist a nces of all 5 SNe via the 
method described in IVinko et all (|201lh . We used two sets 
of velocities for each SNe: i) from the absorption minimum 
of the Fe A5169 line, ii) v mo d e i determined in Sec. 13.41 

The result ed distances are in Table [3] The correction 
factors of lDoBI were applied for all SNe. Usually the photo- 
metric data were interpolated to the epochs of the velocities. 
However, for SN 2004dj Eq.[2]and Eq.|3]were used to extrap- 



olate the velocity data to the photometric epochs, because of 
the low number of spectra taken before day +50, i.e. during 
the expansion of the photosphere. 



6.1 SN 1999em 

Data on SN 1999em in NGC 1637 were used for distance de - 
termination with EPM several times. lHamuv et al.l (|200ll ) 
used cross-correlation velocities (with the model spectra 
of E96 as the template set) and the correction factors of 
E96 to obtain the distance of 7.8 ± 0.5 Mpc. Using ab- 
sorption minima velo cities and the same correction factors, 
iLeonard et al. (|2002l) determined the distance as 8.2 ± 0.6 
Mpc. while lElmhamdi et al.l (120031) o btain ed 7.8 ± 0.3 Mpc. 
On the other hand, ILeonard et al. I (|2003h determined the 
distance of NGC 1637 using Cepheids as 11.7 ± 1 Mpc, 
which is significan tly higher. Using the SEAM method, 
iBaron et aiT(|2004T ) got 12.5 ± 1.8 Mpc. With the velocities 
of th eir CMFGEN models an d and the correction factors of 
D05. lDessart fe Hillierl <|2006h derived 11.5+ 1 Mpc , and with 
a similar approach to that of Baron et al.l ll2004h they ob- 
tained 12.2 ± 2 Mpc. Recently 



Jones et alj (|2009h estimated 



the photospheric velocity from vap (3S3]Eq.[T}, and derived 
9.3 ± 0.5 Mpc by using the correction factors of E96, and 
13.9 ±1.4 Mpc from the correction factors of D05. 

We have repeated the EPM analysis using D05 cor- 
rection factors and our u mo <2ei velocities. This resulted in 
12.5 + 1.4 Mpc, which is in good agr eement with the c e pheid - 
and SEAM-distances and that of iDessart fe Hillierl (|2006T ) 
using the velocities determined from their CMFGEN mod- 
els. Instead, applying the VFe velocities (extrapolating for 
the first few points using Eq[5J), the distance became lower, 
9.7+0.8 Mpc. The disagreement between these two distances 
is roughly the same as that due to the application of different 
correction factors (see above). The distance obtained from 
adopting the v mo d e i velocities is in much better agreement 
with the independent Cepheid-based distance to the host 
galaxy. It suggests that the application of the proper veloc- 
ity data is important to obtain more realistic and bias-free 
distances from Type II-P SNe. 
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6.2 SN 2004dj 

There are many published distances for the host galaxy of 
SN 2004dj (NGC 2403), but they show large scatter, be- 
ing between 2.88 and 6.43 Mpc, according to the NEE0 
database. 

In the case of this SN the velocity curve was extrap- 
olated using Eq. [2] and Eq. p to t he epochs of t h e pho - 
tometry of 1 Vinko et~aH (|2006T ) and iTsvetkov et~all (2008). 
The distances from the two velocity curves agree very well 
(D = 3.6 ± 0.6 and 3.7 ± 0.8 Mpc, respectively, Table [3]), 
and they are also i n very good agreement with the result of 
IVinko et alJ l|2006r i. 



6.3 SN 2004et 

Similarly, the distances of NGC 6946, the host galaxy of 
SN 2004et, show large scatter being in the range of 4.7 
an d 7.2 Mpc (NED). D ~ 4.7 Mpc was derived recently 
bv lPoznanski et al l |2009h using SCM. 

Our result supports this shorter value. The EPM- 
analysis with both v mo dei and VFe (after extrapolation to 
the early phases by Eq. [2]) gave 4.8 Mpc (Table [3]). 



Table 3. The EPM distances of the 5 SNe using different veloc- 
ities. 



6.4 SN 2005cs 

In the case of SN 2005cs the application of v mo del resulted in 
a considerably longer distance (8.6 ± 0.2 Mpc) than the one 
using VFe (7.5 ± 0.2 Mpc). This longer distance is in good 
agreement with a recent study by Vinko et al. (2011), in 
preparation, who determined the distance of M51 via EPM 
by combining the data of SNe 2011dh and 2005cs, and ob- 
tained 8.4 ± 0.7 Mpc. The reason for the slight difference 
between their result and ours with v mo dei (although both 
papers used t he same photometry , velocities and method) 
is that, unlike IVinko et aU (|201ll ). we did not fix the mo- 
ment of explosion in EPM. Instead, we also optimized that 
parameter to keep consistency with the analysis of all the 
other objects in this paper. 



6.5 SN 2006bp 

For SN 2006bp and its host galaxy, NGC 3953, the distances 
are between 15.7 and 21.0 Mpc (NED). Both of our results 
fit into this wide range, but with the usage of v mo dei we 
obtained slightly longer values (20.7 ± 1.8 Mpc) than with 
V Fe (18.6 ± 1.5 Mpc), t he latter being closer to the distance 
of iDessart et~aH l|2008| ), i.e. 17.1 and 17.5 Mpc from SEAM 
and SCM, respectively. 



7 CONCLUSIONS 

In this paper we investigated three methods for estimat- 
ing photospheric velocities of Type II-P SNe. We focused on 
building model spectra with SYNOW, and compared the result- 
ing v m odei velocities with those obtained by cross-correlation 
or simply measuring absorption minima of P Cygni features. 
Based on a sample of 81 spectra from 5 SNe, we showed 



SN D (Mpc) 

with v mode i with V F e 



Photometry 1 



1999em 

2004dj 

2004et 

2005cs 

2006bp 



12.5 (1.4) 
3.6 (0.6) 
4.8 (0.4) 
8.6 (0.2) 
20.7 (1.8) 



9.7 (0.8) 

3.7 (0.8) 

4.8 (0.6) 
7.5 (0.2) 
18.6 (1.5) 



1,2,3 
4,5 
6,7 



1 Source of ph otome try: (1) 
iLeonard et~aI1 i2002T), (3) 
(41 I Vinko et al.l (1200^1. (5) 



(6; 

Pa 



6) Maguirc ct al. 




Hamuv et aU jgOOlf). (2 ) 



Elmhamdi ct al. 



Tsvetkov et al. 

I 



Sahu et alJ d2006h~ (8) 



2003), 
2008), 



http://ned.ipac.caltech.edu/ 



Pastorello et all ( 2009). f9) IDessart et al.l ( 2008) 



that SYNOW provides very similar photospheric velocities to 
those derived by more sophisticated modeling codes, but in 
a faster, less computation-intensive way. This approach may 
be more extensively applicable, yet it preserves the advan- 
tages of using physically consistent model spectra to esti- 
mate parameters of SNe non-interactively, and without re- 
lying mostly on eye-ball estimates and human decisions. 

We illustrated that the cross-correlation- and absorp- 
tion minimum velocities, i.e. those determined by more con- 
ventional methods, suffer from phase-dependent systematic 
deviations from the model velocit ies. This has already bee n 
known from previous studies (e.g. IDessart fc HillierlkoOSbl ). 
but we have extended the phase coverage of the modeled 
spectra, and revealed that such deviations become stronger 
below v p hot ~ 3000 km s _ , i.e. after day +60. At these late 
phases VFe may overestimate v mo dei by 30 - 50 % depending 
on the atmospheric properties of the particular SN. 

Based on these results, we verified and updated the 
relations between the photospheric velocities and the ones 
estimated from the Doppler-shifts of the absorption min- 
ima of individual spectral lines. It was found that while the 
VFe I v mo dei ratio appears to be nearly the same for all SNe 
studied here, it is not true for the velocities from the H/3 line. 
We have derived a power-law relation to estimate v mo dei 
from VFe and/or Vnp, but due to the possibility of SN- 
dependent systematic deviations, we recommend the com- 
putation of parametrized models, whenever possible. 

Using the model velocities, we re-determined the dis- 
tances of the 5 SNe via EPM, and compared them with the 
ones calculated by using vp e - The distances obtained from 
Vmodei are similar or slightly higher than those with vf b . For 
SN 1999em, which is the most thoroughly studied object in 
our sample, we were able to show that by using the model 
velocities the derived distance is more consistent with the 
Cepheid-based distance to the host galaxy. Although such a 
comparison was not possible for the other SNe due to the 
lack of reliable Cepheid distances, this result underlines the 
importance of the velocity measurement method in SN dis- 
tance studies. 

Despite its numerous advantages, EPM also suffers from 
caveats. One of them is the need for many photometric data 
and contemporaneous velocities (i.e. spectra) covering most 
of the plateau phase. This is hardly achievable for most SNe. 
A possible solution may be a careful interpolation between 
the measured data points. Previously, the weakest link was 
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the poorly resolved velocity curve, thus, mainly the light 
curves were int erpolated to the moments of velocity mea- 
surements (e.g. lHamuvl l2002). Based on our results in Sec. 
[4] the interpolation of the velocity curve to the epochs of 
photometric data via Eq.Q may also be a possibility, re- 
sulting in a better sampled dataset for EPM. We intend to 
demonstrate the application of this approach for new SNe 
in a future paper (Takats et al., in preparation). 
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Figure Al. Early spectra of SN 2004et 



the observed spectra of SN 1999em as templates (v cc #l) and 
the CMFGEN models (u cc#2 ). 



APPENDIX C: PARAMETERS OF THE 
BEST-FITTING SYNOW MODELS 

In Table IC1I we present the parameters of the best-fitting 
SYNOW models, including: r re f of the main atoms/ions, the 
power-law exponent of the optical depth function (n), the 
photospheric temperature (Ttb) and photospheric velocity 
(vphot) (see £13.41 for details). 



Table Al. Journal of spectroscopic observations of SN 2004et 



Date 


JD-2,450,000 


Phase(d) 


Airmass 


Observer 1 


2004-10-03 


3281.6 


+11 


1.2 


JT, TK 


2004-10-04 


3282.8 


+12 


2.0 


SM, JT 


2004-10-05 


3283.5 


+13 


1.1 


HD, TK 


2004-10-06 


3284.9 


+14 


2.3 


HD, JT 


2004-10-07 


3285.5 


+15 


1.1 


HD, JT 


2004-10-08 


3286.9 


+16 


1.7 


JT, JG 



1 Observers: JT: J. Thomson, TK: T. Koktay, SM: S 
Mochnacki, HD: H. DeBond, JG: J. Grunhut 



Wang, X., Yang, Y., Zhang, T., Ma, J., Zhou, X., Li, W., 

Lou, Y.-Q., Li, Z., 2005, ApJ 626, L89 
Zhang, T., Wang, X., Li, W., Zhou, X., Ma, J., Jiang, Z., 

Chen, J., 2006, AJ 131, 2245 
Zwitter, T., Munari, U., Moretti, S., 2004, IAU Circ. 

8413, 1 



APPENDIX A: EARLY SPECTRA OF SN 2004et 

Soon after its discovery (|Zwitter et alj [2004), six low- 
resolution optical spectra of SN 2004et were taken at the 
David Dunlap Observatory, Canada, with the Cassegrain 
spectrograph mounted on the 74" telescope. The spectra 
covered the 4000 - 8000 A regime with a resolution of 
~ 800 at 6000 A ( see the details on the data reduction in 
IVinko et all 120061 ). Due to the fixed North-South slit direc- 
tion, the spectra could not be taken at the parallactic angle, 
thus, the slope of the continuum in the blue is affected by 
differential refraction. Table lATI contains the journal of these 
observations, and the spectra are plotted in Fig. I All 



APPENDIX B: SNE VELOCITIES MEASURED 
WITH DIFFERENT TECHNIQUES 

In Table IB1I we present the velocities obtained with SYNOW 
{Vmodei, see t|3.4|) . along with those measured from the ab- 
sorption minima of H/3 and Fen A5169 (fir/3 and i>F e ), as 
well as those obtained with cross-correlation method using 
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Table Bl. The velocities (in km s -1 ) obtained with different techniques: the photospheric velocities from the best-fitting SYNOW models 
(.v mo del)y the Doppler-shift of the absorption minima of H/3 and Fell A5169 features (vjjp and Vp e , respectively) and the cross-correlation 
velocities from templates consisting of the observed spectra of SN 1999em (u cc #i) and the CMFGEN models (v C c#2)- The phase (in days) 
were calculated relative to the date in Table [T] 



SN 1999cm SN 2004et 



phase 




VH/3 






f cc #2 


phase 




VH/3 


VFe 






4.79 


11050 (300) 


11332 






10341 (454) 


11.10 


9700 


(450) 


11072 




9951 (153) 


8781 (1316) 


5.84 


10900 


(350) 


11101 






9858 (888) 


12.30 


8900 


(250) 


8878 




8605 (86) 


8470 (681) 


6.84 


9950 


(250) 


10141 






9112 (878) 


13.00 


9100 


(400) 


9386 




9724 (117) 


8786 (1099) 


7.64 


8900 


(200) 


9148 






8383 (832) 


14.40 


9200 


(400) 


9375 




7617 (71) 


8570 (482) 


8.67 


8850 


(150) 


8697 


- 




8544 (491) 


15.00 


8800 


(100) 


8894 


- 


8635 (64) 


8545 (524) 


12.84 


8550 


(500) 


8406 


7236 




8513 (497) 


16.40 


8400 


(400) 


9282 




9110 (73) 


9060 (545) 


14.14 


7350 


(500) 


7913 


6829 




7575 (289) 


24.60 


7300 


(550) 


8018 


6416 


6624 (49) 




14.67 


8150 


(350) 


7992 


7213 




7708 (190) 


30.60 


6050 


(300) 


7487 


5535 


5779 (56) 




15.14 


7650 


(150) 


7802 


6793 




7380 (187) 


35.50 


5500 


(200) 


7091 


4861 


5112 (43) 


5739 (350) 


19.67 


6750 


(200) 


7246 


6019 




6703 (242) 


38.60 


5100 


(400) 


6869 


4695 


4981 (45) 


5513 (337) 


24.66 


6150 


(300) 


6393 


5146 




5779 (126) 


40.70 


4600 


(600) 


6842 


4644 


4943 (41) 


5448 (256) 


24.84 


6150 


(225) 


6450 


5172 




5866 (194) 


50.50 


3750 


(300) 


5859 


4040 


4236 (49) 


4252 (367) 


27.84 


5000 


(200) 


6060 


4773 




5446 (155) 


54.00 


4050 


(200) 


5468 


3715 


3727 (23) 


4352 (244) 


29.84 


4650 


(250) 


5816 


4693 




5221 (201) 


55.60 


3600 


(450) 


5396 


3686 


3867 (37) 


3895 (408) 


33.84 


4300 


(100) 


5251 


4312 




4394 (215) 


55.76 


3900 


(250) 


5578 


3840 


4048 (37) 


4113 (228) 


37.84 


3900 


(100) 


4885 


4083 




4521 (159) 


63.50 


3200 


(400) 


5018 


3539 


3656 (43) 


3678 (349) 


41.04 


3600 


(200) 


4164 


3421 




3691 (115) 


73.60 


2600 


(450) 


4674 


3292 


3327 (45) 


3793 (217) 


47.84 


3200 


(400) 


3822 


3528 




3765 (163) 


79.73 


2600 


(100) 


4807 


3133 


3192 (31) 


3504 (246) 


50.74 


3400 


(150) 


3699 


3435 




3661 (131) 


84.79 


2400 


(200) 


4725 


3147 


3172 (51) 


3490 (257) 


51.76 


2800 


(350) 




3183 




3147 (152) 


86.50 


2250 


(200) 


4561 


3060 


2876 (64) 


3496 (148) 


66.76 


2300 


(100) 


2645 


2869 




3097 (146) 


98.60 


1900 


(450) 


4377 


2740 


2779 (22) 


3268 (131) 


79.84 


1800 


(400) 


1928 


2313 




1634 (100) 


103.7 


1800 


(150) 


4412 


2779 


2708 (33) 


3060 (201) 








SN 2004dj 












SN 2005cs 






phase 


"^model 


■"Hp 






f cc #2 


phase 


^model 


VH/3 








47.89 


3350 


(300) 


4569 


3183 


3154 (40) 


3471 (83) 


3.44 


7100 (200) 


7275 




17105 (1876) 


30045 (365) 


50.59 


2750 


(150) 


4089 


2920 


2943 (63) 


3401 (172) 


4.41 


6900 (150) 


6848 


6814 


7703 (109) 


7092 (705) 


52.89 


2900 


(150) 


4122 


2949 


3075 (75) 


3243 (115) 


5.39 


6500 (250) 


6487 


6418 


7594 (151) 


6945 (586) 


55.89 


2650 


(100) 


3849 


2897 


2967 (72) 


3124 (122) 


8.43 


5900 (300) 


6050 


5481 


5745 (143) 


5775 (71) 


64.86 


2500 


(150) 


3317 


2625 


2672 (45) 


2887 (114) 


8.84 


5950 (100) 




5222 


5422 (209) 


4337 (193) 


65.85 


2550 


(100) 


3282 


2573 


2678 (47) 


2649 (106) 


14.36 


4150 (300) 


5052 


4167 


4060 (52) 


4119 (80) 


67.86 


2400 


(150) 


3468 


2645 


2700 (46) 


2689 (118) 


14.44 


4600 (100) 


5099 


4213 


3908 (43) 


4591 (170) 


83.89 


2200 


(200) 


2967 


2350 


2384 (50) 


2387 (114) 


17.35 


4150 (50) 


4532 


3763 


3616 (49) 


3757 (129) 


94.67 


1850 


(250) 


3010 


2101 


2222 (67) 


2454 (145) 


18.42 


4150 (100) 


4345 


3788 


3566 (58) 


3742 (110) 


95.87 


1850 


(200) 


3178 


2186 


2302 (71) 


2535 (161) 


22.45 


3550 (300) 


4069 


3049 


3232 (51) 


3765 (315) 


97.87 


1700 


(150) 


3132 


2182 


2270 (78) 


2504 (156) 


34.44 


1900 (50) 


2222 


2318 


2254 (19) 


2505 (126) 


99.86 


1700 


(150) 


3118 


2075 


2208 (86) 


2433 (153) 


36.40 


1800 (100) 


2003 


2142 


2189 (38) 


2202 (74) 
















44.40 


1350 (50) 




1756 


1907 (23) 


2036 (182) 
















61.40 


1050 (50) 




1311 


1459 (24) 


1611 (274) 


SN 2006bp 


phase 


^model 


VHf3 


VFe 


Vcc#l 


V C c#2 


phase 


V-model 




VFe 


Vcc#l 


'"cc#2 


5.30 


12050 (400) 


11251 




11537 (138) 


11131 (571) 


24.26 


7050 (200) 


8278 


6315 


6573 (45) 


7317 (329) 


7.10 


10750 (450) 


11097 




11407 (103) 


10748 (875) 


32.26 


6150 (200) 


6996 


5199 


5548 (28) 


6308 (330) 


9.10 


9600 (200) 


10526 




10975(70) 


10242 (657) 


41.21 


5050 (125) 


6078 


4445 


4698 (46) 


5281 (181) 


11.11 


9550 (450) 


10265 




8717 (69) 


9941 (521) 


56.19 


4050 (200) 


5562 


3952 


4022 (40) 


4602 (180) 


15.10 


8750 (250) 


9735 


8013 


8373 (67) 


9176 (325) 


72.04 


3850 (150) 


5313 


3697 






20.28 


8100 (100) 


9104 


7267 


7506 (79) 


8147 (184) 
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Table CI. The SYNQW parameters of the best-fitting models of SN 1999em. The phases were calculated from the assumed moments 
of explosions listed in Table [T] 



SN 1999em 



phase 


JD-2, 451,000 










're 


/ 










n 


Tbb 


'"phot 


(days) 


(days) 


Hi 


He I 


Nai 


Si I 


Sin 


Can 


Sen 


Tin 


Fen 


Ban 




(kK) 


(km s — ) 


4.79° 


481.79 


2.80 


0.25 


















3.0 


14.2 


11050 


5.84 6 


482.84 


3.50 


0.20 


















3.0 


12.0 


10900 


6.84 6 


483.84 


4.90 


0.40 


















3.0 


11.0 


9950 


7.64 a 


484.64 


6.30 


0.35 


















3.0 


9.5 


8900 


8.67™ 


485.67 


7.30 


0.20 


















3.5 


13.6 


8850 


12. 84 6 


489.84 


15.80 


0.10 


















5.5 


10.0 


8550 


14.14 6 


491.14 


21.10 




0.20 












0.30 




5.0 


11.0 


7350 


14.67" 


491.67 


26.15 


0.05 


0.10 






2.20 






0.80 


0.35 


6.5 


11.5 


8150 


15.14 a 


492.14 


20.20 




0.10 




0.25 








0.70 




4.5 


10.4 


7650 


19.67 a 


496.67 


42.05 




0.10 






16.0 


0.05 


0.30 


0.95 


0.40 


5.0 


9.0 


6750 


24.66 a 


501.66 


113.50 




0.55 






131.9 


0.25 


2.10 


1.85 




8.0 


8.3 


6150 


24.84 6 


501.84 


81.00 




0.35 






17.6 


0.15 


1.10 


1.60 


0.40 


7.0 


8.3 


6150 


27.84 b 


504.84 


73.95 




0.40 








0.30 


1.80 


2.25 


0.20 


6.0 


8.2 


5000 


29.84 b 


506.84 


65.00 




0.45 








0.15 


1.65 


4.10 




6.0 


7.0 


4650 


33.84 b 


510.84 


42.20 




1.10 


0.01 




95.0 


0.05 


1.20 


2.20 




4.5 


5.2 


4300 


37.84 6 


514.84 


39.60 




1.00 


0.01 


0.10 




0.15 


1.55 


2.25 




4.5 


5.5 


3900 


41.04 6 


518.04 


57.20 


0.20 


1.70 


0.01 


0.35 


458.0 


0.20 


4.05 


4.25 


0.20 


7.0 


6.8 


3600 


47.84 i > 


524.84 


23.00 




1.85 


0.01 


0.35 


200.0 


0.45 


2.30 


2.20 


1.00 


3.0 


4.5 


3200 


50.74 a 


527.74 


30.80 




2.05 


0.01 


0.35 


90.0 


0.65 


3.75 


3.35 




4.5 


5.3 


3400 


51.76" 


528.76 


36.60 




1.95 


0.02 




212.0 


0.20 


0.50 


2.65 


0.25 


3.0 


6.5 


2800 


66.76 a 


543.76 


16.80 




5.65 


0.03 


0.40 




1.25 


7.45 


6.45 




3.5 


5.5 


2300 


79.84 b 


556.84 


71.50 




13.95 


0.02 


0.65 


969.0 


2.35 


17.95 


13.90 


4.70 


3.5 


6.0 


1800 


SN 2004dj 


phase 


JD-2,450,000 












/ 










n 




v phot 


(days) 


(days) 


Hi 


He I 


Nai 


Si I 


Sin 


Call 


Sen 


Tin 


Fen 


Ban 




(kK) 


(km s — 1 ) 


47.89 


3234.89 


81.0 




0.70 




0.15 




0.30 


2.85 


2.55 


0.85 


6.0 


7.65 


3350 


50.59 


3237.59 


88.0 




1.10 


0.01 






0.70 


6.15 


6.00 


1.00 


5.5 


6.0 


2750 


52.89 


3239.89 


68.0 




1.05 


0.01 


0.40 




0.40 


3.75 


3.05 


1.00 


5.5 


8.0 


2900 


55.89 


3242.89 


51.0 




1.15 


0.02 


0.20 




0.30 


2.75 


3.70 


0.80 


4.5 


6.2 


2650 


64.86 


3251.86 


72.0 




1.95 


0.01 


0.25 




0.75 


4.25 


4.95 


0.50 


5.5 


8.0 


2500 


65.85 


3252.85 


63.0 




2.05 


0.01 


0.20 




0.50 


4.75 


5.15 


2.00 


5.5 


6.0 


2550 


67.86 


3254.86 


83.0 




2.25 


0.02 


0.30 




0.70 


5.45 


8.55 


0.05 


5.5 


6.8 


2400 


83.89 


3270.89 


85.0 




7.10 


0.01 


0.30 




1.00 


6.30 


7.55 


1.60 


5.5 


8.9 


2200 


94.67 


3281.67 


145.0 




18.45 


0.01 


0.10 




1.00 


7.15 


8.80 


3.55 


5.5 


9.7 


1850 


95.87 


3282.87 


109.0 




24.75 


0.02 


0.30 




1.05 


9.55 


10.20 


2.45 


6.0 


9.4 


1850 


97.87 


3284.87 


185.0 




34.30 


0.06 


0.40 




1.40 


11.70 


14.20 


3.00 


5.0 


9.5 


1700 


99.86 


3286.86 


120.0 




25.0 








0.50 


9.00 


10.0 


0.40 


4.5 


8.5 


1700 



S ource of spectra : 

a |Leonard et al.l &200± 
b iHamuv et alj iboOll ) 
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Table CI - continued for SN 2004et. 



SN 2004et 



phase 
(days) 


JD-2,451,000 

(days) 


LI T 

11 1 


He I 


IN a I 


Ol T 

hi I 


bl II 


/ 

La II 


IT 

OC II 


Ti II 


Fen 


Ba II 


n 


Tbb 
(kK) 


Vphot 

(km s -1 ) 


11.10 a 


281.60 


3.75 


0.20 




















4.5 


80.0 


9700 


12.30 a 


282.80 


4.05 


0.20 




















3.0 


19.0 


8900 


13.00 a 


283.50 


6.40 


0.20 




















5.0 


52.0 


9100 


14.40° 


284.90 


8.00 




0.30 


















4.0 


9.5 


9200 


15.00 a 


285.50 


7.25 


0.15 




















5.0 


43.0 


8800 


16.40 a 


286.90 


14.45 


0.25 




















4.0 


9.0 


8400 


24.60 6 


295.10 


10.0 




0.10 






4.00 









85 


0.20 


3.0 


11.0 


7300 


30. 60 6 


301.10 


21.0 




0.10 




0.40 


74.4 




0.15 


1 


20 




4.5 


7.9 


6050 


35.50 6 


306.00 


41.0 


0.10 


0.05 






97.0 




0.85 


2 


15 


0.05 


5.5 


7.4 


5500 


38.60 b 


309.10 


47.0 




0.25 


0.01 




98.0 


0.15 


1.45 


2 


80 


0.05 


4.0 


7.2 


5100 


40.70 6 


311.20 


36.0 




0.15 




0.05 


95.8 


0.25 


0.75 


2 


05 


0.30 


3.5 


6.8 


4600 


50.50 6 


321.00 


40.0 




0.65 






324.0 


0.25 


1.25 


2 


35 


0.70 


3.0 


5.4 


3750 


54.00 c 


324.50 


81.7 




0.90 


0.01 


0.75 


32.7 


0.55 


4.60 


3 


55 




6.0 


7.6 


4050 


55. 60 6 


326.10 


45.0 




0.95 


0.01 




45.0 


0.45 


3.05 


5 


05 


0.45 


3.5 


6.8 


3600 


55.76 c 


326.26 


58.8 




1.20 


0.02 




25.0 


0.20 


2.90 


2 


50 




4.0 


7.3 


3900 


63.50 6 


334.00 


31.0 




1.15 






81.0 


0.25 


2.0 


3 


60 




2.5 


5.0 


3200 


73.60 6 


344.10 


79.0 




2.30 






100.0 


0.60 


2.75 


5 


05 


1.10 


3.0 


4.8 


2600 


79.73 c 


350.23 


79.9 




4.90 






30.0 


1.20 


11.9 


6 


80 




4.5 


7.1 


2600 


84.79 c 


355.29 


142.0 




7.35 


0.01 


0.50 


10.0 


0.65 


9.15 


5 


15 




3.5 


6.2 


2400 


86. 50 6 


357.00 


223.0 




3.15 




0.10 


515.0 


0.70 


7.90 


13.1 




4.0 


5.2 


2250 


98. 60 6 


369.10 


72.25 




6.05 






50.0 


0.80 


3.10 


4.70 


1.10 


2.5 


4.5 


1900 


103.70 c 


374.20 


219.0 




32.05 






80.0 


2.35 


37.7 


25.3 


5.50 


4.5 


5.5 


1800 



Source of spectra: 

a This paper . 

6 I Sahu et ail feoq6h 

c iMaguire et al. I ll2010bh 



Table CI - continued for SN 2005cs. (The second column is JD-2,453,000). 



SN 2005cs 



phase 
(days) 


JD 

(days) 


III 


He I 


Nil 


Nai 


Mgi 


Si i 


7~ref 

Sin 


Can 


Sen 


Ti II 


Fen 


Ball 


n 


Tbb 
(kK) 


'Vphot 

(km s _1 ) 


3.44 


552.44 


3.90 


0.30 


0.10 




















4.0 


15.0 


7100 


4.41 


553.41 


7.00 


0.40 


0.10 








0.10 












4.5 


15.0 


6900 


5.39 


554.39 


10.6 


0.30 


0.10 




















5.5 


13.6 


6500 


8.43 


557.43 


15.3 


0.20 










0.50 










0.35 


5.0 


10.6 


5900 


8.84 


557.84 


9.00 


0.30 










0.60 










0.45 


9.5 


9.8 


5950 


14.36 


563.36 


81.0 


1.50 




0.40 






0.60 


30.0 


0.10 


0.55 


4.25 


0.30 


8.0 


10.2 


4150 


14.44 


563.44 


49.0 


1.00 




0.25 


0.35 




0.50 


27.0 


0.20 


0.90 


1.30 


0.35 


7.0 


9.8 


4600 


17.35 


566.35 


124.0 






0.60 


0.70 




0.40 


80.0 


0.05 


1.70 


2.05 




8.0 


8.1 


4150 


18.42 


567.42 


168.0 






0.90 


0.50 




0.70 


5.00 


0.50 


1.40 


1.50 




8.5 


9.6 


4150 


22.45 


571.45 


36.0 






0.30 


1.00 




0.40 




0.70 


1.20 


1.05 


0.90 


8.0 


7.5 


3550 


34.44 


583.44 


52.0 






5.10 


0.40 




0.30 


800.0 


2.30 


6.90 


9.0 


4.0 


5.0 


6.5 


1900 


36.40 


585.40 


35.0 






4.75 


0.80 




0.30 


100.0 


1.85 


8.85 


10.3 


1.95 


5.0 


6.2 


1800 


44.40 


593.40 


23.0 






23.9 




0.01 


0.35 


200.0 


4.75 


25.45 


24.1 


7.85 


5.0 


6.2 


1350 


61.40 


610.40 


15.0 






119.9 




0.04 


0.25 


800.0 


12.35 


73.50 


54.4 


63.40 


5.5 


5.7 


1050 
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Table CI - continued for SN 2006bp 



SN 2006bp 



phase 
(days) 


JD-2,451,000 
(days) 


Hi 


He I 


Nil 


Nal 


Sin 


T ref 

Can 


Sc II 


Tin 


Fen 


Ban 


n 


(kK) 


v phot 
(km s -1 ) 


5.30 


840.30 


2.15 


0.15 


















4.0 


10.2 


12050 


7.10 


842.10 


3.20 


0.15 


















3.5 


12.0 


10750 


9.10 


844.10 


6.05 


0.15 


















4.5 


11.7 


9600 


11.11 


846.11 


10.35 


0.20 


















5.5 


12.7 


9550 


15.10 


850.10 


15.65 


0.05 


0.04 


0.20 


0.80 








0.35 




4.5 


10.8 


8750 


20.28 


855.28 


26.70 




0.03 


0.20 


0.40 


23.0 


0.08 


0.20 


1.30 




4.5 


9.5 


8100 


24.26 


859.26 


20.00 


0.10 


0.01 


0.10 


0.10 


25.0 


0.10 


0.80 


1.50 


0.50 


5.5 


8.5 


7050 


32.26 


867.26 


59.70 


1.50 




0.55 




129.0 


0.30 


1.65 


1.75 


0.20 


6.0 


8.0 


6150 


41.21 


876.21 


118.6 


2.00 




1.25 




326.0 


0.65 


3.50 


2.60 


0.85 


7.0 


7.2 


5050 


56.19 


891.19 


145.1 


1.00 


4.65 






300.0 


0.20 


6.40 


5.35 


0.20 


6.0 


7.2 


4050 


72.04 


907.14 


150.1 






18.4 




165.0 


1.90 


10.85 


6.75 


1.45 


5.5 


7.0 


3850 



